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STREAMUNE CURVATXJRE DESIGN PROCEDURE FOR 
SUBSONIC AND TRANSONIC DUCTS 

By Raymond L. Barger 
Langley Research Center 

SUMMARY 

A procedure for designing ducts for subsonic and transonic speeds is described. 
Examples discussed are a wind-tunnel contraction cone, a supersonic nozzle, and a dif- 
fuser. A listing of the computer program is included. 

INTRODUCTION 

When the design parameters (length, maximum diameter, and velocity variation) of 
a duct permit a very gradual area variation, the one -dimensional approximation (ref. 1) 
can be used for tiie design. Such is not the case, however, for many applications, such as 
wind-tunnel contraction sections, nozzles, and engine inlet diffusers for which there are 
severe restraints on allowable dimensions. For ducts of this type the radial variation of 
the flow parameters is often not negligible, and so the one- dimensional approximation does 
not suffice. 

If a duct is intended to be used for low-speed applications, it can be designed by a 
generalization (see ref. 2 and its references) of an incompressible flow method originated 
by Tsien (ref, 3). The usual procedure for the compressible case is to utilize a series 
expansion in powers of some parameter (refs. 4 to 7), but such solutions rapidly increase 
in complexity with increasing number of terms, and therefore must be truncated after a 
small number of terms. For example, the solution of reference 5, which utilizes an 
expansion in powers of r2, becomes unmanageably complicated even for the first five 
terms. 

The problems associated with series expansion solutions are avoided by the method 
of reference 8, which can be used as a design procedure. That method is a relatively 
complicated direct numerical attack on the problem, involvii^ a continuous transforma- 
tion of coordinates and a finite-difference solution of the differential equations working 
outward from the center line. 

The present analysis, in comparison, makes use of the streamline curvature equa- 
tions , which represent a natural formulation of the equations for this type of problem 
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inasmuch as any streamline can be taken to be a wall contour. The streamline curvature 
equations have also been used for the design of centrifugal impellers (ref. 9). The pres- 
ent method differs from that of reference 9 in several respects: 

(1) The design velocity distribution is specified on the axis rather than on an inner 

hub profile 

(2) The equations are formulated directly in terms of the streamline curvature 

rather than the radius of curvature, which becomes singular at an inflection 

point of a streamline 

(3) The graphical method of reference 9 is replaced by a mathematical model which 

can be programed for machine calculation 

The streamline curvature equations, as used in the present analysis, represent a 
form of the exact, compressible, inviscid flow equations. The method is therefore appli- 
cable from low subsonic to supersonic speeds, but no treatment of shocks is included. The 
equations, and the computer program in the appendix, have been written for the axisym- 
metric case, but the analysis can be applied to problems of two-dimensional design sim- 
ply by replacing the factor 27rr in the continuity equation with the width of the duct. Then 
a symmetric duct can be designed by specifying a velocity distribution along a center line, 
or one wall can be specified together with a velocity distribution on the wall. 

SYMBOLS 

The International System of Units (SI) is used in this study, 
a streamline curvature, d<^/ds 

C contraction ratio 

c speed of sound 

M Mach number 

m mass flow rate 

n distance normal to streamlines 

q total velocity 

s distance along streamline 
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V 


center-line design velocity 


X,R cylindrical coordinates nondimensionalized with respect to estimated entrance 

diameter 

x,r cylindrical coordinates 

y ratio of specific heats 

p density 

I dr 

— for a curve r = r(x) 


Subscripts: 

c center line 

i ith point on a given streamline 

k kth streamline 

t total conditions 

wall at wall 

A bar over a symbol indicates average value: 


<)k>n 

2 


k+1 


ANALYSIS 
Flow Equations 

When the circumferential component of velocity is neglected, the compressible 
inviscid streamline curvature equations for axisymmetric flow consist of the continuity 
equation in integrated form, 

r*n 

m = 27 t J pqr dn (1) 
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and the velocity equation, 


q = qc 



( 2 ) 


where the density is given by 


,l/(y-l) 


p-p. 


1 - 


Y - 1/q 


and the curvature by 

a. 32 

ds 


(3) 


(4) 


(See ref. 9.) The streamline curvature is the rate of change of flow ai^le with distance 
alor^ the streamline; it is the reciprocal of the radius of curvature. 

Method of Solution 

The calculation is initiated by prescribii^ a velocity distribution along the center 
line. For purposes of calculation the velocity on the kth streamline is determined by the 
approximation 




aAn 


(5) 


The mass flow between streamlines is estimated by 


Am = 2vrpq An 

In order to compute the curvature a^, A<p is determined as the change in ai^le of slope 
of the straight lines joining the ith point to the i - 1 point and to the i + 1 point and 
then A(p is divided by the corresponding As. 

The equations were solved by an iterative numerical technique. First, the shape of 
a stream tube having a small maximum diameter is estimated from the one-dimensional 
flow relations usii^ the assumed center-line velocity distribution. The normals to this 
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streamline are determined, as well as its slope and curvature distributions. These func- 
tions permit the velocity and density distributions to be computed on this streamline by 
means of equations (5) and (3). Next, the mass flow is computed through each normal 
and compared with the required mass flow for this stream tube. The streamline coordi- 
nates are then adjusted slightly in the direction that will reduce the error in mass flow, 
the error is computed, and the coordinates are readjusted, and so on until the desired 
mass flow is obtained within a reasonable tolerance. This procedure is carried out at 
each point on the streamline. Then the curvature of this streamline is computed and used 
to obtain a new estimate for a to be used in equation (5). With this new estimate the 
entire procedure is repeated to obtain a better estimate for the streamline. The velocity 
along this streamline is then taken as input data to compute the next streamline, which is 
first estimated from the one -dimensional equations, in the same manner as for the pre- 
vious streamline. 

The successive streamlines are computed similarly until a specified mass flow is 
obtained, unless the solution becomes unstable before the mass flow condition is met. In 
such case the last stable streamline can be taken as the wall contour, unless the solution 
displays an unacceptable pressure gradient along the wall. 


Stability of Convergence 

One drawback to the use of the stream curvature equations is the tendency of the 
solutions to become unstable relatively easily. This tendency is a direct result of the 
use of the streamline curvature. A small chaise in r can result in a relatively large 
charge in the curvature. Since the curvature is in the exponent in equation (2), the veloc- 
ity is sensitive to the curvature and thus to small changes in r. 

This sensitivity can be counteracted to a considerable extent by the manner in which 
the calculation is performed. In the present method, a smooth distribution of curvature 
is used in the initial estimate for each streamline. Then when the actual curvature is 
computed, it is not used directly but is averaged with the curvature of the previous 
streamline for use in equation (5), This procedure tends to stabilize the convergence and 
has been used in the calculations of this paper. 

There is an alternate procedure which is mathematically more stable, but it 
requires somewhat more computing time. If many streamlines are computed, the follow- 
ii^ approximation will apply in equation (2): 

^ ^n 

r a dn = ^ a^(An)^ - I 
" 1 1 
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Thus the coordinates of each streamline can be computed from the distributions along 
the previous streamline, and it is not necessary to determine the coordinates of a given 
streamline from its own curvature by iteration. 

Design Velocity Distributions 

Some of the existing design analyses (e.g., refs. 3 and 5) prescribe the design axial 
velocity distribution in terms of a particular analytic expression. Where high-order 
derivatives are required, this analytic expression must be selected so that these deriv- 
atives do not become unmanageably complicated. 

The present method, on the other hand, is not restricted to a particular type of 
design velocity distribution inasmuch as the input data for the calculation can be in the 
form of values specified at discrete x-locations on the axis. This system facilitates 
"experimenting" with various designs by makii^ small local variations in the design dis- 
tribution. However, the input data should be reasonably smooth in order to obtain wall 
shapes without excessive curvature, as well as to Insure stable convergence of the cal- 
culation. For some designs requiring large values of the stream function and rapid flow 
changes, the design velocity must be extremely smooth. This requirement can be met by 
inputtir^ a smooth second derivative and integrating it twice. 

The selection of a design velocity distribution need not be simply a cut-and-try 
process. The one -dimensional equations are often useful to provide an estimate of the 
flow relations. Moreover, in those cases where both the upstream and downstream flows 
are uniform, the one-dimensional equations give correct relationships between the flow 
values in these two regions, even though the flow between these regions is nonuniform 
(provided, of course, that no separation takes place). 

For example, suppose that a wind-tunnel contraction cone is to be designed for an 
exit Mach number of 1 and a contraction ratio of C. Then the uniform entering velocity 
would be determined by the equation 



(See ref. 10, eq. (82).) 

Since in the enterii^ flow V/c^ is generally much less than 1, the nonlinear term 
can be neglected in an approximation, and so, for y - 1.4, v/c^ = 125/216C. The veloc- 
ity at the exit would be given by (ref. 10, eq. (68)). 
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Examples 


Three sample cases have been computed: a wind-tunnel contraction section, a noz- 
zle, and a diffuser. The contraction section is shown in figure 1 together with several 
streamlines. The most critical part of such a duct is that near the exit, where the flow 
becomes parallel and, correspondii^ly, the wall curvature becomes zero. A wall which 
is otherwise smooth but has a rapid decrease in curvature incurs an adverse pressure 
gradient at the wall, with the possibility of flow separation and the certainty of a severe 
radial velocity gradient. For the example shown, the curvature variation is gradual, so 
that the flow is fairly uniform radially, as indicated. 

The other critical region for a contraction section is near the wall at the large end, 
where the slow-moving flow enters a region of negative curvature. This adverse gradient 
may separate the flow locally. However, for the sake of obtaining a short contraction 
section, some separation might be tolerated in this region with the expectation that the 
subsequent favorable acceleration section will reduce flow nonuniformities. 

A nozzle design is shown in figure 2. Observe that where the wall curvature is 
slight, the variation of Mach number along a streamline normal is slight, although the 
duct may be sharply contracting. On the other hand, the radial variation is relatively 
large in the throat, where the cuiwature is greater. The supersonic Mach number dis- 
tribution at the wall is smooth because of the smooth variation of wall curvature. This 
result may be compared with that obtained with a cone — circular-arc — cone nozzle, which 
displays a drop in wall Mach number at the curvature discontinuity in the expanding sec- 
tion of the nozzle (see ref, 11, fig. 1). 

Figure 3 shows a design for a subsonic diffuser. This case is more sensitive than 
that of the contraction section because the negative curvature at the wall occurs where 
the flow is already decelerating. Since this critical region is near the terminal end of 
the diffuser, the results of any flow separation will be seen at the end of the diffuser. 

To illustrate the influence of the curvature on the solution, the wall Mach number 
distribution for the nozzle is plotted in figure 4, as calculated from both one -dimensional 
and streamline curvature equations. The differences are larger in the vicinity of the 
throat, where the Mach number is relatively sensitive to flow geometry. The maximum 
difference in Mach number is about 0.06. 

CONCLUDING REMARKS 

A computerized procedure that makes use of the inviscid streamline curvature 
equations has been developed for designing compressible-flow ducts from an assumed 
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center -line velocity distribution. Three examples are given: a wind-tunnel contraction 
section, a nozzle, and a diffuser. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 4, 1973. 
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APPENDIX 


COMPUTER PROGRAM 

The program language is FORTRAN 4. The input, by the NAMELIST format, con- 
sists of the following quantities: 

X an array of axial positions for which velocities are specified 

Q an array of the design velocity values 

RMX the desired radius at X(l) 

ROI Pt 

CT c^ 

KMX number of streamlines to be computed 

N number of X and Q values input 

Before beginning the calculation the program prints the design center-line velocity 
distribution with the headings 

J point index 

X axial coordinate 

Q velocity 

After computing each streamline the program prints the results for that streamline 

with the following headii^s: 

PT, point index 

ITER. number of iterations required to satisfy continuity at Jth point 

X axial coordinate 
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APPENDIX - Continued 


R 

radial coordinate 

SLOPE 

streamline slope 

CURV. 

streamline curvature 

VEL. 

velocity along streamline 

P/PT. 

local pressure divided by total pressure 

M 

local Mach number 


Computing times for the three sample cases ranged from 2 to 5 seconds. 


PROGRAM MfiKADllCi INPUT, 0UTPUT«TA»P‘ I=IMPIJT> 

MlMpNSIOl^ X(50> ’O(SO) *SL<S0) *RrS0> «XJ<S0) ,P.J(G0> (SO > » SK (SO) * 

1^0 (SC) .GK| (SO) *01 (SO) 

NAMrLTST/lNPT/X ,0*RMX,R0T «M«CT.KMX 

C f'EbJ) TNOUT ANO RPfNT VPLOCTTY niSTPlRUTlOM, MONnlMENSlONALlZE VfL, 

1 RPADd ,TNPT)$tF(E0F,1)?7»? 

? Rr = 3.14lSQ ‘BKrl 
mKtmT 40 

i,r, format ( IHl *RX1HJ, HX1 hX,1 I XIHQ/) 

PRIMT 41 , (J,X (J) ,0(0) ,.J=1 ,N) 

00 6 J=].M 1iO( J) ='M J)/CT $SKI(J)=,0 1>0I ( J) =0 ( J) %SK(J)=.Q 

41 P OPMAT ( T I0»?n 2*G) 
fs P I ( . J ) = . 0 

K OP KTm STRFAmL tNP; , COMPUTE PAPAMFTFRS AMO TNTTtAL FSTIMaTFS fOR 
COORD TNaTeS AMD SLOPES 

7 kITsELOAT (K)«RMX/EL0AT (KMX) SAI=PT»(PT I*RTT-PI (D »HI (1 ) ) 

‘^M = ROI*Al*On)»(l,-.2*Qn)*0(I) ) <s DO R J=1 ,M S0RS=0 ( J) <»0 { J) 
PO(J)=POI<M1.-.'^-»ORS)**2,S % A = Q (1 )*AI/(RO( n «0(J) ) 

4 K ( ji=SORT (A/PJ+RI (J)*RI ( )> ) *5Ln)=,0 SOO <5 J=P,N 

C COMPUTF normals to kTm STRFAMLINf 

R SL (J)=(R(J)-R(J-1) )/(X(J)-X(J-l) >S no lA J=?*N 
IF(J.F0.N)G0 to 13 

iF (SL (J) .r,T.,g, AMO, SL(J*1 ) ,GT, 0)00 TO 10 
1E(SL(J) .LT..O.*UO,SL(J+1>.LT.O)GO TO 11 
IF (SL ( J) *SL(J*1) .LF,.0)OO TO 13 

10 SM=SL(J) SR1=R(J-1) SX1=X(J-1) SPL=R(J)SXL=X(J) SGO TO 1? 

11 SM = SL(J*l) SRl = R(J) $Xl = X(J)i;XL=xrj*-l)SRL=P(J+l) 

1 ? IF (SM.MF. .0) XJ ( J) =SM/ ( SM«SM,l . )*(RI ( J)-Rl*SM<»Xl*X ( J) /SM) 

ie(sm.fo,, j)Xj(J)=xi 

RJ( J) =P1 *SM<^(XJ( J)-X1 ) SGO TO 14 

13 RJ(J)=R(J)S XJ(J)=X(J) 

14 CONTINUE SXJ(1)=,L SRJ(1)=RII SSK(1)=,0 
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APPENDIX - Concluded 


c COMPUTF VFLOCTTY, DFNSTTY, AMO MASS PLOW 

LL = 1 

15 00 17 J=e.N Sl.= l 

16 IF (LL.FO.l . AND. K. GT. 1 > SK (. J) =SK T ( J)*FLOAT(k: ) /float (K-1> 

OM-sORT ( (XJ(J)“X(Jn*»*?*(PJ(J)-RT( J) 1 »•»?> SSKM-fSK (J) +SKT (J1 )*.5 
U< J)=OI ( J) *FXP(SKM*0N1 

0M= (01 < J) ♦Ol J) SROM=( **3,5 

RM=,5* <PJ ( J> +RI ( JM SFMM=?,*PT*RM*nN*ROM*OM5 FR= (fM-FMM ) /FM 
IF(ARS«FR) •l.T,, jOuOOODGO TO 17 

ADJUST COoROlNAfFS OF KTh STRFAMLTNF ALONG NORMAL TO SATISFY 

continuity 

RJ(J)=RT (J) ♦(1.*F‘^*.'5)»(RJ(J)-RI ( J) ) 

IF (ABS (0< J5**0( J-l) ) ,LT,.OaO0OOl 1 R J ( J) =R J ( J-U 
IF (J,EO.?>RJ( J-i)-RJ( J) 

XJ(J)=X (J)+<l,*tR*,S)*(XJ(J>-X(J) > SL=L+1 SIF(L,GT.RR1G0 TO ?7 
GO TO 16 

17 CUNTTNUF 

PRINT 44* < J*X < J) «RJ( J) «SL( J> ♦SK(.n ♦ J= 1 .N) 

44 F0PMAT(i7,4F1?.S) 

kF.COMPUTF CURVATURE FOR AOJUSTFD COOROINATFS ANo RFRFAT 

calculation for kth streamline 

LL=LL*1 $IF(LL.GT.?)G0 TO ?5 
00 IB J=2.N 

IP SL( J) = <Rj(J)-RJ<J-in/<X.HJ)~XJ(J-l) 1 
N1=N-1 500 ?0 J=?*N1 

IF ( (R.J(J)-RJ(J-n )(J>1 .GF..01GO TO *R 

RJ ( J) =,5* (RvH J*1 » +RJ(J-1>> 

SL ( = (RJ{ J*1>-RJ' J-i) )/(X J( J*11-XJ( J-U > 

19 IF (J,LT,?>5L(J)=,v 
?0 CUNTTNUF SSK(1)=,CI SSK(N)=,0 
00 33 J=?*N1 

^>=,5* (SORT ( (Rj( J) -RJ( J-1 ) 1 **?♦ (X J( J)-X J ( J-1 U *SORT ( (RJ( J*1 ) 

I-RJ(J))P*3+<XJ(J*l)-XJ(J)>**?>) 

OFL=* ATAN(5L { J»l) >-AT aN(SL( J) 1 
SK ( J)=(nEl /S*SK( J> )*,5 

33 CONTINUE SOO 34 J=3.Nl S TE ( SK ( J> *5K ( J- U ,LT, . 0 . aNH . 5K ( Jl *SK ( J* 1 > . 
1LT,.0)SK(J)=,S*(SK(J-1) *-SK(J*l) ) 

34 CONTIMUF sSK(l)=,y SSK <3>=SK (3> *,5 SGO TO IS 
!'I0NDTMFMSI0NALI'?F AMO PRINT RFSULTS FOR kTH STRrAMLTNF* SHIFT 
INDEX* AND PROCEED TO NFXT STREAMLINE 

3F PRINT 43 

43 format < 1H1 ♦SXSHITFR, *SX3HPT. *6X1HX. 1 3X1 HR.RXSHSI OPF , lOKSHCUOV, * 

1 I 0X4MVFL . ,9X4hP/PT « 1 3X IhM/> 

00 76 J=3*N 

0T(J1=0(J) SX(J»=:XJ(J) SSKI ( J)=SK( J) 

^0P=1 .-.3»Q( J) *0( J1 TXM=0(J>/SORT(ROP> 5POPT=poP**3,5 

SR=RJ< J>/(3.*RMX) $ SX=X J(J>/r3.*RMX) 

43 f 0 RMaT(3TP*7F13.6> 

print 43*L*J»SX*SR.SL(J» »SK(J) ,0( I) ,POPT, XM 
IF ( K.FO.KMX) 0 ( J) =. 0 

PI(J1=RJ(J) SOKl)=PJ(l) SK=K*1 STF (K.LF.KMXIGO TO 7 5G0 TO I 
?7 STOP send 
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Figure 1.- Wind-tunnel contraction section with streamlines and design velocity distribution. 
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Figure 2.- Nozzle with streamlines and design velocity distributions. 
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Figure 3.- Diffuser with streamlines and design velocity distribution. 
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Figure 4.- Mach number distribution along wall of nozzle shown in figure 2, as obtained by streamline 

curvature calculation and by one -dimensional theory. 
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